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We review the proposal to resum the physics of hard thermal loops in the ther- 
modynamics of the quark-gluon plasma through nonperturbative expressions for 
entropy and density obtained from a <t>-derivable two-loop approximation. A com- 
parison with the recently solved large- Nf limit of hot QCD is performed, and some 
updates, in particular on quark number susceptibilities, are made. 

1. Introduction 

Even at temperatures several orders of magnitude higher than Aqcd the 
perturbative series for the thermodynamic potentials of hot QCD 1 ' 2 ' 3 does 
not appear to converge and thus seems to be devoid of predictive power. 
While the first correction at order g 2 gives a reasonable estimate for the 
results obtained in lattice gauge theory a few times above the transition 
temperature, everything breaks down as soon as collective phenomena such 
as Debye screening come into the play and produce formally higher-order 
contributions suppressed by only single powers of g (see Fig. 1). 

Since similar difficulties have been observed in simple scalar models 6 ' 7 , 
the reason for this failure does not have to do so much with the fact that 
QCD has a nonperturbative sector even in its deconfined phase which limits 
the number of computable perturbative coefficients. It rather seems that 
screening effects should better not be treated in a strictly perturbative way. 
A first encouraging attempt for improving the situation was put forward 
by Karsch et al. 8 , who proposed to keep a screening mass unexpanded at 
any given order of the loop expansion and to fix this mass by a stationarity 
principle. Technically, this corresponds to rewriting the Lagrangian as 

£ = £ - imV + A„t + fm^ 2 . (1) 
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Figure 1. The poor convergence of thermal perturbation theory in pure glue QCD. The 
various grey bands bounded by differently dashed lines show the perturbative results to 
order g 2 , g 3 , g 4 , and g 5 with MS renormalization point p, varied between -kT and A-kT. 
The thick dark-grey line shows the lattice results from the Bielefeld group 4 ; the lighter 
one behind that of a more recent lattice calculation using an RG-improved action from 
the CP-PACS collaboration 5 . 

This optimization of thermal perturbation theory was successfully ap- 
plied to <f) 4 theory to three-loop order 9 and extended to QCD by Andersen 
et al. ' . There they proposed replacing the simple mass term by the 
gauge-invariant hard-thermal-loop (HTL) action 12 ' 13 . This method, which 
they termed HTL perturbation theory (HTLPT), is in principle a system- 
atic and manifestly gauge invariant scheme. From a physical point of view, 
it has, however, the somewhat unsatisfactory feature that HTL's are used 
uniformly for soft and hard momenta although the HTL effective action is 
accurate only for soft momenta (and soft virtuality) <C T. In the following 
we shall aim at a resummation that captures more closely the physics of 
the collective excitations of the quark-gluon plasma. 

2. ^-derivable approximations and HTL resummation 

Our proposal 14 ' 15 ' 16 is to implement HTL resummations in the spirit of the 
so-called ^-derivable approximations 17 of the 2PI skeleton expansion. In 
the latter, the thermodynamic potential is expressed in terms of dressed 
propagators (D for bosons, S for fermions) according to 

Q[D,S] = ±TTr logD" 1 - ^TTtUD 

-TTr logS^ 1 + T Tr + T$[D,S] (2) 

where & [D, S] is the sum of 2-particle-irreducible "skeleton" diagrams. The 
self-energies n = D^ 1 — D^ 1 and S = S^ 1 — S^ 1 , where D and 5*0 are bare 
propagators, are themselves functionals of the full propagators, determined 
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by the stationarity property 

5Q[D, S}/6D = = 5Q[D, S}/6S. (3) 

according to 

5$[D,S]/5D=±I1, 6$[D,S\/6S = S. (4) 

The ^-derivable two-loop approximation consists of keeping only the 
two- loop skeleton diagrams, which leads to a dressed one- loop approxi- 
mation for the self-energies (4). In a gauge theory this introduces gauge 
dependences (which are however parametrically suppressed 18 ), but we shall 
construct further approximations which are manifestly gauge independent. 

A self-consistent two-loop approximation for fl has a remarkable conse- 
quence for the first derivatives of the thermodynamic potential, the entropy 
and the number densities: 



dT 



9, T ^ 

Because of the stationarity property (3), one can ignore the T and /x depen- 
dences implicit in the spectral densities of the full propagators, and differ- 
entiate exclusively the statistical distribution functions n and / in (2). Now 
the derivative of the two-loop functional T$[£>, S] at fixed spectral densities 
of the propagators D and S turns out to cancel part of the terms Im (HD) 
and Im (SS 1 ) in (2) leading to the remarkably simple formulae 19,15,16 

5 - tr /(0^[ Iml °^"- ImnRe ^ 

- 2tT J (ty^ [imlogfi^-ImEReg], (6) 
AA=-2tr| I g I ^[lmlog^ 1 -Im S Re^]. (7) 

Through this formulae, all interactions below order g 4 are summarized 
by spectral data only, which shows that entropy and density are the pre- 
ferred quantities for a quasiparticle description. 

The leading-order interaction terms cx g 2 arise from hard loop momenta 
involving only the light-cone projection of the self-energies, e.g. in pure glue 
QCD 

-Tre(u>)8(u> 2 -k 2 ) 
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where Ily is the transverse component of the gluon self-energy, which on 
the light-cone is accurately (to order g 2 ) given by its HTL value lir(fc, k) = 
™oo = even though k is no longer soft 20 ' 21 . 

The more critical g 3 ( "plasmon" ) term in the thermodynamic potentials, 
on the other hand, arises in an unusual manner: while in the pressure it is 
determined by the soft momentum regime of the one-loop contribution, in 
the above expressions for the entropy, there are two distinct origins. One 
part comes from order-g corrections to Hr(k, k) at hard momenta k ~ T 
and are to be interpreted as a correction to the entropy of hard gluons. 
Only a fraction (< 1/4) arises as the entropy of soft gluons in the HTL 
approximation 14 . The required correction for the hard excitations is 

and this can be calculated by standard 22 HTL perturbation theory 3 -. 

These corrections to the asymptotic thermal masses are, in contrast to 
the latter, nontrivial functions of the momentum, and can be evaluated 
only numerically However, as far as the generation of the plasmon term is 
concerned, these functions contribute in a certain averaged form which can 
be calculated analytically, 

- 2 f dkkn'(k)ReSU T (uj = k) 1 2 _, 

*»~ - - — — = -^ 9NTmD - (l0) 

This result pertains only to the hard excitations; corrections to the vari- 
ous thermal masses of soft excitations are known to differ substantially from 
(10). For instance, the relative correction to the gluonic plasma frequency 23 
at k = 0, Snipt /m 2 ; , is only about a third of Sm^/m 2 ^; the NLO correc- 
tion to the nonabelian Debye mass on the other hand is even positive for 
small coupling and moreover logarithmically enhanced 24 . 

For an estimate of the effects of a proper incorporation of the next- 
to-leading order corrections we have therefore proposed to include the lat- 
ter only for hard excitations and to define our next-to-leading approxima- 
tion (for gluons) through S NLA = S H TL\ soit + | hard , where m 2 ^ in- 
cludes (10) and where Shtl refers to evaluating the ^-derivable 2-loop 
entropy in the HTL approximation exactly. To separate soft (k <~ Tod) 
and hard (k ~ 2irT) momentum scales, we introduce the intermediate scale 
A = \Z2-KTmDCK and consider a variation of ca = \ . ■ . 2 as part of our 
theoretical error estimate. 



a Bccause the external momentum is hard, no HTL vertices are needed, and only one 
propagator has to be dressed by HTL. 
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Another crucial issue concerns the definition of the corrected asymptotic 
mass fhoo- For the range of coupling constants of interest (g > 1), the 
correction |5to^J is greater than the LO value m^, leading to tachyonic 
masses if included in a strictly perturbative manner. 

However, this problem is not at all specific to QCD. In the simple g 2 (p 4 
model, one-loop resummed perturbation theory gives 

m 2 =g 2 T\l- *g) (11) 

which also turns tachyonic for g > 1. On the other hand, the solution of 
the one-loop gap equation is a monotonic function in g, and it turns out 
that the first two terms in a (m/T)-expansion of this gap equation, 

m 2 = g 2 T 2 _ Ig^Tm, (12) 

7T 

which is perturbatively equivalent to (11), has a solution that is very close 
to that of the full gap equation (for MS renormalization scale p, ~ 2ttT) 16 . 

In QCD, where the non-local gap equations are too complicated to be 
attacked directly we adopted (12) as a model to include Smi^, which is 
needed for the completion of the plasmon effect. b This leads to 

™ 2 oo = U N + \N f )g^ -^-g^NTm^. (13) 

Up to a single integration constant the resulting entropy expres- 
sion determines the thermodynamic pressure. Choosing this (strictly 
nonperturbative 14 ) input such that e.g. P{T C ) « 0, where T c is taken from 
the lattice, this ambiguity in fact becomes quickly negligible for larger T/T c , 
because the contribution of the bag constant thus introduced drops like T~ 4 
in the normalized pressure P/Po, where Po is the ideal-gas limit. 

The main uncertainty rather comes from the choice of the renormaliza- 
tion point p,. In the following we always consider varying fi by a factor of 2 
around a central value of 27rT and determine the strong coupling constant 
from the 2-loop renormalization group equation. 

The result of this procedure is displayed in Fig. 2 and compared with 
lattice data from the Bielefeld 4 and PC-PACS 5 groups, and also with the 
recent two-loop calculation in HTLPT 11 . This shows a clear improvement 
compared to the perturbative result to order g 3 , and remarkable agreement 
with lattice data for T > 3T C . 



b Numerically, this differs only slightly from the Padc approximants that we employed in 
our first publications 14,15 . 
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Figure 2. The pressure in pure-glue QCD in the next-to-lcading approximation devel- 
oped in Ref. 16 (BIR) compared to perturbation theory through order g 3 and to the 
recent 2-loop HTLPT result 11 (ABPS), all with p varied between ttT and 4nT. The 
gray bands represent (continuum extrapolated) lattice results 4 ' 5 . 



Compared to the HTLPT calculation, an important difference of our 
approach is the separate treatment of hard and soft contributions, but the 
HTLPT result also has a large ^-contribution 11 with opposite sign from 
that obtained in 3-loop perturbation theory. 



3. Inclusion of fermions and the large- Nf limit 

When fermions are included (Nf ^ 0), part of the plasmon effect in the 
pressure (and all of the plasmon effect in the fermion density) is contributed 
by next-to-leading order corrections to the asymptotic thermal mass of 
the fermions, whose leading-order (HTL) value is — g 2 CfT 2 /4, with 
Cf = 2N/N g . These can be calculated in standard HTL perturbation 
theory according to 

^SMl(k) = <5S + (^ = k)= Re( , £X + ./*\ , (14) 

which is again a function of k that can be evaluated only numerically. 
However in the plasmon effect it enters only in the analytically calculable 
averaged form 

- SM 2 jdkk.f>(k)Re2k6Z + (u = k) 1 

SM °° ~ Jdkkf(k) -2^ 9 CfTm ° ■ (15) 

In our previous work 16 we have incorporated this correction in complete 
analogy to the gluonic asymptotic mass (13), that is, we have adopted a 
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gap equation quadratic in M x that is perturbatively equivalent to (15). 

However, the recent work on the laxge-Nf limit of QCD by Moore 25 
has triggered us to reconsider this procedure, because it turns out that a 
quadratic gap equation for the fermions does not comply with the large- iV/ 
limit. 

In the large- Nf limit, Nf — > oo, g 2 — > such that g 2 s = g 2 Nf/2 <~ 1, 
the quadratic gap equation for the gluons has the correct behaviour that 
ml^g^Tye + Oil/Nf). 

Fermion self-energies are suppressed by a factor 1 /Nf , but we still need 
to consider them in our expressions for entropy and density because there 
are Nf fermions which together produce a Nf contribution to the fermionic 
entropy. This precisely equals —NNfM^T/6, where represents the 
average appearing in (15). The latter involves the correction term calcu- 
lated in (15), but without further (rainbow-like) corrections on the internal 
fermion line. The fermionic "gap equation" thus has to remain linear in 
. An equation for Ml with the correct behaviour in the large- Nf limit 
is given by 

Ml = g 2 C f T 2 /A - -^g 2 C s Tm 00 . (16) 

There is then no negative feedback from the fermion mass itself, it only 
inherits higher-order terms from the solution to rhoo (when Nf is finite). 

In Fig. 3 we compare the 2-loop <I>-derivable result in our approximations 
in the limit of large Nf with the exact result c . The latter has the curious 
behaviour of being nonmonotonic as g 2 s is increased, and this behaviour 
is in fact qualitatively reproduced in our approach with (16). However, 
when g 2 e gl 7.4, Ml becomes negative (dashed lines in Fig. 3). This is 
not a problem for the ^-derivable expressions at order N®, but it means 
that for large but still finite Nf the fermionic quasiparticlcs (at least in our 
approximation) cease to exist. 

However, in our applications to real QCD the revised fermionic gap 
equation (16) does not have the problem of giving rise to tachyonic masses 
even close to the transition temperature. If we therefore compare the results 
of our approximations in the large Nf limit with the exact one only in the 
region where Ml remains positive, the outcome is in fact encouraging: the 
agreement below the point where Ml vanishes is remarkably good even 
though the coupling is no longer small and rho/T = g c tf./V3 <~ 1. 



c The result published in Ref. 25 has recently been found to be in error; the corrected 
exact result can be found in Ref. 26 . 
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By a curious coincidence, for N = 3 and Nf = 3 the revised gap equation 
(16) together with (13) has exactly the same solutions as the uncoupled 
quadratic gap equations we used previously. Only for Nf > 3 there is 
at all a coupling where the fermionic mass ceases to grow monotonicly 
with g; for Nf < 3 this never happens. Because of this coincidence, the 
numerical changes in our previous results are almost completely negligible. 
As an example, Fig. 4 updates the results published previously 16 for the 
entropy with flavour numbers Nf — 0, 2, 3 in comparison with the estimated 
continuum extrapolation of Ref. 27 . Only the NLA result for Nf = 2 
changes at all from the gray dash-dotted line to the one slightly below it. 



P/P (large-TVf limit) 
1.1 



S/S (N f = 0, 2, 3) 




Figure 3. The exact result 26 for the 
pressure in the limit of large Nf com- 
pared with the ^-derivable 2-loop result 
in the HTL approximation (full lines) 
and in the next-to-leading approximation 
(full lines ending in dashed lines), for 
p, = T and 4irT. The gray lines denote 
the next-to-leading approximation with 
quadratic fermionic gap equation consid- 
ered in Ref. 16 , but which we argue needs 
to be replaced by Eq. (16). 



Figure 4. The entropy in the HTL ap- 
proximation (full lines, Nf = 0, 3, 2 from 
bottom to top) and in the next-to-leading 
approximation (dash-dotted lines, Nf = 
0, 2, 3 from bottom to top) with fi=2irT 
and compared with an estimated contin- 
uum extrapolation 27 of the lattice result 
for N f =2. Only the NLA N f =2 result 
changes by switching to the fermionic gap 
equation (16) from the gray dash-dotted 
line to the black one just below. 



More interestingly, with the new fermionic gap equation (16) we can 
narrow down somewhat our predictions 28 for the quark number suscepti- 
bilities. Previously we have determined our estimated theoretical errors for 
the latter by combining the results obtained from a quadratic gap equation 
with the results from a Pade approximant. To be compatible with the large- 
Nf limit, we can now restrict to (16) and produce somewhat narrower error 
bands (which still are dominated by the p, dependence). Fig. 5 shows our 
prediction for Nf = 2 in comparison with available lattice data 29 . Fig. 6 
shows the results for quenched QCD together with recent results for two 
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different continuum extrapolations 30 (both are higher than the previous 
lattice results 31 ). 

It would be interesting to compare our predictions with HTLPT to two- 
loop order; the one-loop HTL-rcsummed results from charge correlators 32 
suffer from severe overcounting while missing out the plasmon term 33 so 
that one should not perform a comparison yet. d 

x/xo 





Figure 5. Comparison of our results for 
X/XO in masslcss Nt=2 QCD with the 
lattice results of Ref . 29 (no continuum ex- 
trapolation). Full lines refer to the HTL 
approximation, dash-dotted lines to NLA. 
(Gray dash-dotted denotes our previous 
NLA estimate 28 .) 



Figure 6. Comparison of our updated 
results for x/xo hi the formal limit 
Nj=0 with the previous lattice results 
for quenched QCD of Ref. 31 (gray 
data points) and two recent continuum 
extrapolations 30 (black data points). 



4. Outlook 

The HTL-based quasiparticle description of the thermodynamics of hot 
QCD that we have developed can be straightforwardly extended to finite 
chemical potential. First steps in this direction that go beyond the simple 
quasiparticle models of Ref. 36 are indeed encouraging 37 . Further refine- 
ments, in particular a full inclusion of the momentum dependence of the 
next-to-leading order asymptotic masses, is work in progress. 
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